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We construct an effective model for the chiral field and the Polyakov loop in which we can investigate the 
interplay between the approximate chiral symmetry restoration and the deconfinement of color in a thermal 
SU(3) gauge theory with three flavors of massive quarks. The phenomenological couplings between these two 
sectors can then be related to the recent lattice data on the renormalized Polyakov loop and the chiral condensate 
close to the critical region. 
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I. INTRODUCTION 

Various theoretical developments accompanied by several 
numerical results from lattice simulations have advanced our 
understanding of QCD thermodynamics JTHj]. Nevertheless, 
there are still some important qualitative questions that remain 
to be answered. For instance, as one increases the tempera- 
ture above the critical one, which mechanism drives the QCD 
phase transition for a given quark mass? Is the QCD tran- 
sition more deconfining or chiral symmetry restoring? Al- 
though lattice QCD might provide definitive responses soon, 
some current findings are still controversial. In particular, the 
numerical value of the critical temperature for the QCD phase 
transition yD, the uncertainty on the location and even the 
possible absence of a critical point in the temperature-baryon 
density plane [4], and the newly suggested possible difference 
between the critical temperatures for the chiral and and the de- 
confinement transitions |5|]. Therefore, it is useful to address 
these issues from a phenomenological point of view, resorting 
to effective field theory models. 

The pure gauge sector of QCD with N colors, correspond- 
ing to infinitely heavy quarks, is well under control in lattice 
simulations, with a precise prediction for the deconfinement 
critical temperature, and a good understanding of its thermo- 
dynamic behavior [6]. This sector has a global Zn symmetry 
associated with the center of the gauge group SU(N). The 
Polyakov loop i, charged under Z#, serves as order param- 
eter for the deconfining transition. It is real for N — 2, and 
the Z2 symmetry breaking deconfinement transition is of sec- 
ond order [7]- For N = 3, £ is complex and transforms as 



„2ni 



Accordingly, the expectation value (£) van- 



ishes at low temperatures, when Z3 is unbroken, and (£) ^ 
above the deconfinement critical temperature T c / ~ 270 MeV. 
This transition is verified to be weakly first order |J. 

Another theoretically well-studied sector of QCD is the 
limit of zero and small quark masses, i.e. the chiral limit. 
In this regime Zn is always broken, but chiral symmetry is re- 
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stored above a critical temperature T t . The order parameter 
is the quark chiral condensate a. The transition is believed 
to be of second order for two llfj, ll ill , and of first order for 
three flavors of massless quarks. For finite quark masses chiral 
symmetry is explicitly broken. Increasing m q from zero corre- 
sponds to weakening the first-order phase transition. For some 
quark mass the first-order critical line turns into a second- 
order chiral critical point. For even larger m q the transition 
becomes a crossover. Effective field theories have been used 
to analyze this sector lfl3ll . 

Deconfinement and chiral symmetry restoration are phe- 
nomena of very different nature, corresponding to approxi- 
mate symmetries of the QCD Lagrangian. For realistic quark 
masses both chiral and Zn symmetries are explicitly broken, 
and most likely both transitions are in the crossover domain 
IU4I1 . The phase diagram illustrating the quark mass depen- 
dence of the deconfinining and the chiral symmetry restoring 
transition has been qualitatively investigated in effective the- 
ories in [H d] and on the lattice in d [H [H. In the 
crossover region neither a nor £ can play the role of a true or- 
der parameter. However, lattice simulations show that the sus- 
ceptibilities associated with these two quantities peak at the 
same temperature when quarks are in the fundamental repre- 
sentation of the gauge group [20], suggesting that the tran- 
sitions occur at the same critical temperature. Similar results 
were obtained also in terms of the chemical potential 12111 . Re- 
cent, more refined results confirm the coincidence of the two 
transitions |22|], although results using a different method find 
critical temperatures that differ by ~ 25 MeV |5j], showing 
that this issue is still not settled. 

There exists quite a substantial effort to understand the true 
nature of these transitions u sing different effective field theo- 
ries m m H3, H Hi . We especially consider Ref. 
j27tl . in which the authors provided an explanation, within a 
generalized Ginzburg-Landau theory, to why for m q = chi- 
ral symmetry restoration leads to deconfinement. The analysis 
was based on the following general idea: the behavior of an 
order parameter induces a change in the behavior of non-order 
parameters at the transition 13011 via the presence of a possi- 
ble coupling between the fields, gi£o 2 . In 1I27I1 it was assumed 
that g 1 > 0, which is in line with recent results from the lattice. 

In this paper we extend the analysis of 12711 to the case of 
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nonzero quark masses. Within this model we can then not 
only study the qualitative relationship between deconfinement 
and chiral symmetry restoration including massive quarks, but 
also discuss how to provide a few semi-quantitative predic- 
tions for the behavior of the condensates and critical tempera- 
tures adopting values for the couplings constrained by lattice 
data. 

Lattice QCD results for the temperature dependence of 
the chiral condensate for three degenerate flavors of massive 
quarks were reported in 03111 . and the first results for the be- 
havior of the renormalized Polyakov loop under the same 
lattice conditions were presented in 03211 . Most recent lat- 
tice data, using improved lattice actions, provide the temper- 
ature dependence of the chiral condensates IT331 l34ll and the 
Polyakov loop for 2+1 flavors IT341, l35ll . They indicate that the 
chiral and deconfinement transitions coincide but, contrary to 
the case of the chiral condensate, the flavor and quark mass de- 
pendence of the Polyakov loop is contained entirely in the fla- 
vor and quark mass dependence of the critical temperature. In 
l33tl the strange quark condensate was studied separately from 
the light quark condensate, and it was found that chiral sym- 
metry restoration happens together for the quarks of different 
flavors. Furthermore, it was shown in 03411 that, as expected, 
the strength of the chiral transition decreases while that of the 
deconfining transition increases with increasing quark mass. 

This paper is organized as follows. In Section II we de- 
scribe the effective theory in which we couple the Polyakov 
loop to the chiral sector. In Section III we discuss how to re- 
late our model to the recent findings from the lattice in order 
to constraint the couplings of our model. Here we set up the 
theory. We will report on the quantitative analysis elsewhere 
ll36ll . Section IV contains our summary and outlook. 

II. THE EFFECTIVE THEORY 

For nonzero quark masses both the chiral and the Z# sym- 
metries are explicitly broken, so that neither the chiral conden- 
sate nor the Polyakov loop play the role of a true order param- 
eter any longer. One can, however, still use these quantities to 
monitor the approximate transitions, as is done in lattice QCD 
calculations, where the transition temperatures are determined 
from the peak position of the corresponding susceptibilities. 
Using these quantities we can thus construct an effective field 
theory to study the interplay between the approximate chiral 
symmetry restoration and the deconfinement of color. In what 
follows, we write down the most general renormalizable ef- 
fective Lagrangian that can be built from the chiral field G and 
the Polyakov loop £, following the form introduced in 12711 . 

For the chiral part of the potential we chose the Si/ (3) x 
SU (3) linear sigma model 

V x (4>) = m 2 Tr[* t 4>]+A, 1 (Tr[(I> t 4>]) 2 +A,2Tr[(O t O) 2 ] 

- c(det[«I>]+det[* t ])-Tr[//(* t +#)], (1) 

where 4> = T a § a — T a (a a + i% a ) is a complex 3x3 ma- 
trix containing the scalar o„ and the pseudoscalar K a fields, 
and H = T a h a is the explicit symmetry breaking term. Here 



T a = X a /2 are the nine generators of SU (3) with X a the Gell- 
Mann matrices, A-o = y / 2/31, and h a nine external fields. The 
different symmetry breaking patterns of this theory have been 
thoroughly investigated in [37] and [38]. We are interested in 
making a direct comparison with lattice results both for the 
Nf = 3 degenerate mass flavors, m u = m c i = m s , and for the 
more realistic Nf = 2 + 1 case, with m u = m c i < m s . 

Here we set up the degenerate Nf = 3 case, for which the 
potential (fTJ undergoes some simplifications. This particular 
case corresponds to that named 4a in 13711 . Accordingly, c = 0, 
since the U(1)a anomaly is not taken into account, and ho^Q 
is the only symmetry breaking term, since = h% = 0. In this 
case, the scalar and pseudoscalar matrices are diagonal and 
the masses are given by 

X 

ml = m 2 + 3(A,i + y)a^, 
m% = m 2 + (X l +X 2 )a 2 07 (2) 

7 7 /*\ ^2x9 

< = m- + (Xi + —)ao. 

In the vacuum, the condensate is Go = y/2/ 3f n and ho = 
y273/jim 2 , thus <t> = f n and H = f % m 2 . 

The idea now is to write this special case of the SU (3) po- 
tential in the form 

V x (a) = '^a 2 + ^-Ha. (3) 

This choice is motivated by the fact that the up-down conden- 
sate is aligned in the a direction. Furthermore, without the 
I/(1)a anomaly the difference between the results obtained 
for SI/ (2) and SI/(3) in the sigma model are negligibly small. 
This is clearly illustrated in the melting of the condensates 
shown on Fig. 6 from l38ll . In (fJJ the coupling is given by 
X = X\ + X%/3. The parameters of the model m 2 , Xi and X2 
are determined from the vacuum values of m a , m n and mf . 
In what follows we assume that pions and /o's will not play a 
role, so we can discard them from the outset. 

We adopt a mean field analysis, and thus replace the fields 
by their expectation values. Then, as customary, we choose 
the vacuum to be aligned in the sigma direction. As a result, 
all our relevant equations will depend only on (a) . In the a 
direction the potential has the form shown in Eq. (0). 

The Polyakov loop potential when the Z3 symmetry is ex- 
plicitly broken reads 

Effective field theories using this potential have been intro- 
duced in i39ll and further analyzed in f4Cj[. Here we focus 
only on the real part of the Polyakov loop, since one can al- 
ways choose the expectation value to be real, at least for /j = 0. 
In this case, the contributions coming from the Polyakov loop 
potential simplify to 

V P dt)=8ol+'^e 2 + 8 -±f + 8 -±t*. (5) 
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Here, £ is a scalar field (with dimensions of mass) proportional 
to the Polyakov loop (a dimensionless quantity), and ;«o is its 
bare mass above the transition. The couplings g, that appear 
in each piece of the complete potential can be fixed by lattice 
data as will be shown in the next section. 

The interactions between the chiral field and the Polyakov 
loop which exist when both symmetries are spoiled by the 
presence of massive quarks yield the following contribution 
to the total effective potential 

y mf (£,4>) = ( gl £+g 2 e 2 )TT[^^} 

+ (git 2 +g 2 i)'Tr[H(& + $>)}, (6) 

which simplifies to 

V int (£,a) = gi£o 2 + g 2 £ 2 a 2 + g l £ 2 a + g 2 £a. (7) 

We now expand the fields around their mean values, a = (o) + 
8o and £=(£)+ 8£, and compute the mean-field equations of 
motion from first variations of the complete effective potential 
V — Vpl + V x + Vim, obtaining: 

2X(o) 3 - m 2 a (a) + H-g l (£) 2 - g 2 (£) = (8) 

and 

So + m 2 (£) - g3 (£) 2 - 2g 4 (£) 3 + gl (a) 2 + g 2 (a) = . (9) 

The mass of the o field, m a , is determined as the second 
derivative of the potential with respect to a evaluated at the 
minimum: 

m 2 a = m 2 + 3X(a) 2 + 2 8l (£) + 2g 2 {£) 2 . (10) 
Similarly, the mass of the Polyakov loop is given by 

m\ = m 2 + 2g 3 (£) + 3g 4 (£) 2 + 2g x (a) + 2g 2 (a) 2 . (11) 

^From now on, we neglect higher-order contributions from 
(£), but not from (a). We have to keep higher-order terms in 
the chiral field to be able to recover the correct chiral limit 
discussed in Ref. 12711 . when quarks are massless and H = 0. 
Moreover, couplings which were included in the interaction 
potential of the present analysis due to explicit chiral sym- 
metry breaking must vanish for H = 0. For this reason, we 
require these couplings to be linear in H ~ m q + o(m?). Ac- 
cordingly, we make the following replacements: g\ — » g\H 
and g 2 — > g 2 H. The equations of motion thus become 

21(a) 3 -m 2 a (a)+H-g 2 H(£) =0 (12) 

and 

g + m 2 (£)+g l (o) 2 +g 2 H(o)=0. (13) 
To first order in m q , we have: 

{i) ~ _«o_^( a )2_*2 H ( ff) (14) 
nq rrq 



and 

Notice that in the chiral limit, i.e. for H = 0, we recover the 
results from Ref. l27ll . One can now use ( fT~4-b and ( fT3b to 
extract the couplings from lattice results for the renormalized 
Polyakov loop and the chiral condensate in the presence of 
massive quarks. In the following we discuss how that can be 
implemented in a general fashion rf36tl . 

III. CONNECTING TO LATTICE DATA 

Our aim is to use the above derived equations to describe 
the corresponding lattice data for the chiral condensate j3lll 
and the Polyakov loop 113211 for three degenerate flavors, gen- 
erated under the same lattice conditions. We can apply the 
same procedure for the 2+1 flavor most recent lattice data 
lHH HH HH. in this case we also study separately the light 
and strange quark condensates |[36ll . 

In order to compare our mean-field results to lattice data, 
we use the Gell-Mann-Oakes-Renner (GOR) relation 

w 9 (w) = -/ 7[ m^(a), (16) 

which in the vacuum yields m c/ (ijA|/}o = —f^_m\, and define 
the dimensionless quantity that is measured on the lattice x = 
(W)/(w)o. where we take (w) = -2(225 MeV) 3 . The 
mass of the sigma meson is given by 

ml(T,m q ) = ml(m q )+2'k(o) 2 {T) . (17) 

Furthermore, as customary in a Ginzburg-Landau theory, we 
assume 

ml{T,m q = 0)=aT 2 (l-y) , (18) 

below the transition, where a is a dimensionless constant. 
This is characteristic of a second-order phase transition. For 
nonzero quark mass the transition is a crossover. We then con- 
sider that the temperature and quark mass dependence of the 
sigma mass are well described by the following ansatz 

m 2 a (T,m q ) = m 2 K (m q )+aT 2 (\ - L \ . (19) 

To describe the temperature-dependence of the Polyakov 
loop mass, a fit to pure SU(3) Yang-Mills results from the 
lattice 14111 was performed in [40]. Accordingly: 

m 2 dT)=(^ff ) (20) 

where the temperature-dependent string tension is given 
by S (T) = 1.21^/4 -0.99T 2 ct /(0.41) 2 , with the zero- 
temperature string tension sq — (440 MeV) 2 . This 
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parametrization, however, is not valid in full QCD. Therefore 
we keep the Polyakov loop mass as a free parameter. 

Given the conditions above, we can rewrite ([T5l in a more 
convenient form: 



Co 



g\g2 



mj(T) 



m q x 2 - C\m 2 a (T ,m q ) . 



gogl 
m\{T) 



0, 



(21) 



where we have defined the following constants: C2 = 
-(vjA|/) /2A,/2 ~ 33 MeV, Q = 1 /2kf% ~ 2.9 x 1(T 6 MeV~ 2 
and Co = -(tjty)o/2X/* ~ 3.8 x 1(T 3 MeV- 1 . 
We can also define the dimensionless quantity 



{£) 



(-go/mj(T)) 

gl /„\2 g2 



(-go) 



{-go 



-H(a). 



(22) 



to be compared with lattice results for the renormalized 
Polyakov loop as a function of temperature for different quark 
masses ll32tl . Notice that y has the expected behavior in the 
limit of very high temperatures. In terms of the dimensionless 
chiral ratio x and the quark mass m q : 



1 



£1 



(-So) 



gi 



-go) 



{\pif)om q x. 



(23) 



The functions x = x(T) and y = y(T) for a few values of m q , 
and also the corresponding values of T c , were computed on the 
lattice |[3ll[32ll . By fitting equations (fJTJ and ( |23l , one can ex- 
tract from the data the couplings a, go, g\ and g2, to be used 
in the effective field theory. Our preliminary results indicate 
that the sign of the coupling g\ complies with the expectation 
of l27ll . assuring that the chiral and the deconfinement transi- 
tions coincide. Detailed results with further conclusions are 
deferred to our upcoming publication 13611 . 



massive quarks. We devoted special attention to the behavior 
of the quasi order parameters a and I with temperature, which 
can be connected to lattice data. The latter, on the other hand, 
can be used to provide constraints on the couplings of the ef- 
fective theory. A detailed analysis of the method, as well as 
the extraction of the range of physical values of the couplings 
will be presented elsewhere IBol . 

The effective theory presented in this paper is very simple. 
For instance, if one studies the renormalization of Polyakov 
loops, one is naturally lead to consider effective matrix mod- 
els for the deconfinement transition, which unfolds a rich set 
of possibilities, especially at large N l42ll . In particular, eigen- 
value repulsion from the Vandermonde determinant in the 
measure seems to play a key role l43ll . These studies pointed 
out however, that in the neighborhood of the transition the rel- 
evant quantity is the trace of the Polyakov loop, underlying the 
relevance of our effective theory in this region. The simplicity 
of our model allows for a direct comparison to lattice data in 
a way that provides definite constraints on the couplings, after 
which one can test and predict semi-quantitatively the behav- 
ior of the approximate order parameters in different settings. 

The full determination of the couplings of an effective 
model linking confinement and chiral symmetry breaking us- 
ing reliable lattice data is useful for the study of the phase 
diagram for strong interactions, and can bring understanding 
to some of the open questions considered in the introduction. 
Furthermore, the effective potential derived from this model 
can be used in the study of the real-time dynamics of the QCD 
phase transitions, including the effects from dissipation and 
noise that result from self and mutual interaction of the fields 
related to the chiral condensate and the Polyakov loop 114-411 . 
Results from a real-time study within our effective model will 
be reported in the future [45]. 
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IV. SUMMARY AND OUTLOOK 

In this paper we have sketched a systematic method to con- 
struct a phenomenological generalized Ginzburg-Landau ef- 
fective theory describing simultaneously the processes of chi- 
ral symmetry restoration and deconfinement in the presence of 
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